Natural multiparticle entanglement in a Fermi gas 
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We investigate multipartite entanglement in a non-interacting fermion gas, as a function of fermion 
separation, starting from the many particle fermion density matrix. We prove that all multiparticle 
entanglement can be built only out of two- fermion entanglement. Although from the Pauli exclusion 
principle we would always expect entanglement to decrease with fermion distance, we surprisingly 
find the opposite effect for certain fermion configurations. The von Neumann entropy is found to 
be proportional to the volume for a large number of particles even when they are arbitrarily close 
to each other. We will illustrate our results using different configurations of two, three, and four 
fermions at zero temperature although all our results can be applied to any temperature and any 
number of particles. 
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Introduction. Quantum entanglement plays a crucial 
role in quantum mechanics, and is extensively used in 
quantum information. However, it was only recently 
that researchers have started to investigate entanglement 
in systems containing a large number of particles. This 
is of fundamental importance because entanglement was 
found to be relevant not only in microscopic systems, 
but also on a macroscopic scale P, . Multipartite en- 
tanglement seems to play an important role in condensed 
matter systems, and might be the key ingredient to the 
solution of unresolved physical problems such as high 
temperature superconductivity |3J. In this work we in- 
vestigate multipartite entanglement in a non-interacting 
Fermi gas. Bipartite entanglement in this simple quan- 
tum system has already been shown to be fully charac- 
terized by the exchange integral due to the antisymmetry 
of the wavefunction 0, 0, . 

In this letter we will show that all multiparticle entan- 
glement can be built only from bipartite entanglement. 
This is a significant result because it shows that a com- 
plete description of quantum correlations at all levels is 
possible in a realistic many body system such as a non- 
interacting Fermi gas. We show that the n-particle den- 
sity matrix can be written as a sum of the completely 
mixed state and a mixture of all possible two fermion 
antisymmetrized wave functions. We use entanglement 
witnesses to illustrate that genuine tripartite entangle- 
ment does not exist in this system, in agreement with our 
previous expansion of the density matrix. We then inves- 
tigate bipartite entanglement for three and four fermions 
for different fermion configurations. This entanglement 
is quantified using the negativity Q. Finally we show 
that for large number of fermions the entropy is always 
proportional to this number (which in turn is propor- 
tional to the volume of the system), independently of 
the fermion distance. For a small number of particles 
and small fermion separation the entropy is smaller than 
this number (volume). This clearly establishes the fact 



that entanglement of a non-interacting Fermi gas can be 
treated like any other macroscopic physical quantity and 
that it can be related to other macroscopic observables 
such as the volume of the gas or number density. Any 
mean field theory ignoring entanglement when describ- 
ing macroscopic effects in many body systems is there- 
fore unlikely to be successful even when, remarkably, the 
constituents of the system are non-interacting as in our 
case. 

Density matrix. We consider a many fermion system 
with a fixed number of particles and a density matrix 
p. The elements of the reduced density matrices for 
l,2,3...n particles labeled by pi, p2, p3-..p n respectively 
are given by [8|: 

<*t(iO*(i)> 

(«' t (2')* t (l')*(l)*(2)) 
(«' t (3')* t (2')* t (l , )*(l)*(2)«'(3)) 
(&(n')&((n - l)')-.*(n - l)*(n)) 



(l|pi|l'> 
(12|p 2 |l'2') 
(123|p 3 |l'2'3') 
(l...n\p n \l'...n'} 



where 1 = (ri, o~i), ri is the position vector and o\ =t, \ 
is the spin of the fermion. The average is given by 
(...) = Tr{p...}. For the sake of simplicity all our results 
are illustrated at zero temperature, where p = \4>o)(4>o\j 
with |0 O ) = Ilfe F C I a\ vac ) ec l ua l s the ground state of the 
Fermi system. The ct is the creation operator that cre- 
ates an electron of momentum k and spin a. The Fermi 
momentum is denoted by fcp- and the vacuum state is 
\vac). The '5 are the field operators and obey the usual 
fermion anti-commutation relations , (r^), \& CTl (ri)} 
= S a ' liai S(ri - ri). 

After a somewhat lengthy but straightforward calcu- 
lation we arrive at a form for the density matrix for n 
particles which is particulary useful to investigate entan- 
glement: 



Pn 



(1 



(i) 



2 



where |*~) = l/v2(| TJ-) — I IT)) i s the maximally en- 
tangled singlet state of the pair ij. The sum runs over 
all the pairs ij. The probabilities py are functions of 
the relative distances between all pairs. As an example, 
we write down the density matrix for the two and three 
particle case: 

92 = + |f-)<*- 1 

Pi = (1 -Pl2 -PX3 -P23)^ +Pl2l*i2)(*12l ® \ 
+ Pl3|*r 3 )<*r3l ® 2 +^2 3 |*23><*23l ® ^ 

(2) 

where p = (2 - 2/(r) 2 )/(2 - /(r) 2 ) and /(r) = ji^/a;, 
with the Bessel function j\(x) = (sinx — xcosx)/x 2 and 
x = kpr. The relative distance between the fermion 
pair is denoted by r. The function /(r) is one for 
r = and zero for large r. For three fermions, we 
have three different pairs and for the pair ij : pij — 
+ fijfikfjk)/(-2 + % + fl + fj k -fijfikfsk). The 
function /y is a function of the relative distance between 
fermion i and j only. Note that the probabilities py can 
be calculated for any number of particles. 

Entanglement. The Peres-Horodecki criterion is the 
condition for the existence of entanglement in the two 
particle case 0. In our earlier work, we found this to 
imply that f(r) 2 > 1. This means that two electrons 
are entangled if the relative distance between them is 
smaller then 1.8/fep for T — 6]. Two fermions are 
maximally entangled if they are at the same position. 
This is because of the Pauli exclusion principle. In gen- 
eral, since the overall state must be antisymmetric, if 
the spatial wavefunctions fully overlap and thus are sym- 
metric, then the spins must be antisymmetrised. The 
fermions must, therefore, be in the maximally entangled 
spin singlet state We will show that such two- 

particle entanglement is also the main building block for 
multi-particle entanglement. In order to illustrate this 
behavior we will now consider entanglement in systems 
containing three fermions . 

Tripartite entanglement. From the decomposition of 
the density matrix it is clear that no genuine tripartite 
entanglement exists. We now formally show this using 
the method of entanglement witnesses. These are ob- 
servables which (by our convention) have a positive ex- 
pectation value for all separable states, and a negative 
expectation value for some entangled states, i.e entangle- 
ment exists if Tr{p^H} < 0, where II is the witness 
It has been shown that there are only two different classes 
of tripartite entanglement, which are represented by the 
\GHZ) = 1/V2(| TTT) + I 1U>) and \W 3 ) = 1/V3(| TIT 
) + I ITT) + I TTD) states |l(J. The corresponding wit- 
nesses are defined as: Hghz =1/2— \GHZ)(GHZ\ and 
Uw 3 = 2/3 - |W 3 )(W 3 |. For both of these witnesses, 



because < \fij\ 2 < 1, the trace of p^U cannot be neg- 
ative. This confirms that genuine tripartite entanglement 
does not exist in the ideal Fermi gas. 

Bipartite entanglement . We will now investigate if 
there is entanglement between two groups of fermions. 
One group contains fermion i and the other group the 
fermion pair jk. As a measure of entanglement we use 
the negativity, defined by Nuj k ] = (||/9 3 *||i — l)/2, where 
|| ||i is the trace norm of the partial transpose of the 
reduced density matrix p 3 of fermion i versus the other 
two jk and denote it as Nuj k ] ■ The trace norm can be 
evaluated to be HP3* || 1 = 1+2| Pi I > where the sum goes 
over the negative eigenvalues of the partial transpose. 
There are eight eigenvalues in total and only two of them 
are negative having the same value A. The negativity is, 
therefore, Num = 2|A|. Negativity is a good measure 
of bipartite entanglement because it is monotonic under 
local operations and classical communication and is also 
equal to zero if fermion i is not entangled to the fermion 
pair jk. It reaches its maximal value of 1/2 if fermion i 
is maximally entangled to the fermion pair jk. We will 
now investigate different arrangements of three fermions, 
and investigate the behavior of negativity. 

We first consider three fermions on a straight line. 
The distance between the fermion i and the fermion j 
is fixed. The remaining fermion k moves from the posi- 
tion of fermion i, to the position of fermion j, i.e from 
x = to x m ax — kpfij, where is the relative position 
between fermion i and fermion j. At x — 0, fermion k 
is maximally entangled to fermion i, and can therefore 
not be entangled to fermion j, independent of ry. The 
state of the total system is then p 3 = |^JI)(^'^.| ® \ - As 
fermion k moves away from fermion i we expect the neg- 
ativity N\k t ij] to first drop but and then to increase again 
as fermion k approaches j. This is confirmed in Fig.l. 

We next consider the case where the fermions are lo- 
cated on the edges of an isosceles triangle. Fermions i and 
j form the base of the triangle which is fixed. Fermion 
k is moved away from the midpoint of the base. The en- 
tanglement negativity N^ k ^ and Nuj k ^ for this scenario 
are plotted in Fig. 2. The entanglement N[ k i j] monoton- 
ically decreases as k moves away from ij because the ef- 
fect of antisymmetrization becomes weaker with the dis- 
tance (dashed line in Fig. 2). The entanglement negativ- 
ity Nu jk] (solid line in Fig. 2) initially follows the same 
trend as ^V[fe,ij] for exactly the same reason. Surpris- 
ingly however, the entanglement Nuj k ] , after reaching 
its minimum value, starts to increase and then reaches 
its saturation value. The reason for this is the follow- 
ing. When fermion k is further away from i and j than 
the distance between fermion i and j itself, the effect of 
antisymmetrization between i and j on entanglement is 
larger than the effect of antisymmetrization between i 
and k. If then the distance is further increased, the posi- 
tion of fermion k has a vanishingly small role on entangle- 
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FIG. 1: The negativity ATr 2j i3] is plotted for fermion 2 moving 
from fermion 1 to fermion 3. This corresponds to the text 
values i = l, j — 3, k — 2. The distance between 1 and 3 is 
fixed to Xmax = 5. This is a dimensionless number. 

merit. The three particle density matrix then becomes: 
P3 ~ P% ® 1/2, where P2 is the reduced density matrix 
of the pair ij. Note that the minimum of negativity is 
reached when the fermions are equally distant from each 
other. This is because entanglement is monogamous and 
each particle has to share entanglement equally with the 
other two. This case will now be analyzed in more detail. 

We now consider the case when the fermions are sepa- 
rated by equal distances. For three particles the fermions 
are located on the edges of an equilateral triangle, and 
for the four particle case on the edges of an tetrahedron. 
Entanglement in this case is plotted Fig. 3. We start by 
putting the fermions in a very small volume of radius 
e. Because of the Pauli exclusion principle, only two 
fermions can be in the same location. More then two 
fermions would mean that at least two quantum num- 
bers are the same, which is forbidden. As the distance 
between the fermions increases, Nuj^i] < ^[i,mn] < ^[p.q] 
at all distances, because the entanglement is shared be- 
tween the fermion pairs, and the more fermions are in- 
volved the less entanglement we gain. All of this can be 
generalized to an arbitrary number of fermions. If the 
fermions are all in a small volume of radius e, then the 
state is in an equal mixture of all the singlet states of the 
pairs, and only the second term of Q survives. Higher 
order entanglement does exist, but the Pauli exclusion 
principle forbids maximally entangled states other then 
the |'F - ). This is the reason why we do not have \GHZ) 
or IW3) states in the system. If the fermions are fur- 
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FIG. 2: The negativities ATr 2j i3] = iVr 3il2 ] (solid), JVp. 231 f° r an 
isosceles triangle, for which fermions 2 and 3 form the base 
are plotted. Fermion 1 is moving away from the midpoint of 
the base. This corresponds to the text values i = 2, j = 3 
and k = 1. 

ther away from each other the first term in JTJ becomes 
important and the total state p n becomes even less entan- 
gled. The interplay between the two terms in the density 
matrix is also important when we want to calculate the 
total entropy of the fermions S(n) = —Tr{p n In p n }. We 
study this quantity, because at T = it quantifies the 
amount of entanglement between the measured electrons 
and the remaining unmeasured electrons. 

Entropy. We now investigate the von Neumann en- 
tropy of the Fermi system as a function of fermion dis- 
tance as shown in Fig. 4. At T = 0, the system is in 
a pure state. For the two particle case, the entropy is 
zero for x — because the two fermions are in the pure 
singlet state. It then increases to two as the distance 
increases. This behavior is also observed for three par- 
ticles as well as four. In these cases it does not start at 
zero because we cannot have more then two fermions in 
the same location, but again reaches the value of three 
(four) as the distance increases. For large n, the entropy 
becomes nearly equal to the number of particles even 
for very small distances. For a large particle system at 
low constant density, the reduced density matrix of the 
system is then given by p n — 1/2™ and the entropy is 
S(n) = nln2 which is proportional to the volume of the 
system. For a dense system the density matrix is given 
by: Pn = I ® I/2"- 2 . The entropy of this 

state is also proportional to the number of fermions for 
large n. Only if the number of fermions is small and 
if they are very close, the Pauli exclusion principle pre- 
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FIG. 3: The negativities Mi2]j ^[1,23] (solid) and 
A^[ lj2 34] (dashed) are plotted. This corresponds to text values 
p, q = 1, 2, I, run = 1, 23, and i, jkl = 1, 234. 
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FIG. 4: The yon-Neumann entropy Si (dashed), S3 (solid), S4 
is plotted as a function of fermion distance. This behavior is 
analytically explained in the text. 

vents the entropy from being proportional to the number 
of fermions. This can be explained as follows: If all the 
states are equally likely the entropy is the logarithm of 
the number of all possible configurations of these states. 
Since we have a system of fermions they have to be an- 



tisymmetrized and therefore this number is equal to the 
total number of states 2™, minus the number of symmet- 
ric states (n + 1) . The entropy therefore is: 

S(n) = ln(2" - (n + 1)) (3) 

We can see a complete agreement between this formula 
and the entropy in Fig. 4 for x — 0. It is also clear that 
if the number of fermions is large, then the first term in 
the entropy dominates, giving us the previous result of 
entropy being proportional to the volume of the system. 

Conclusion. We have presented a form of the density 
matrix for n fermions in an ideal fermi gas which is par- 
ticulary useful to investigate entanglement in this system. 
We then showed that no genuine multipartite entangle- 
ment exists, and that all multipartite entanglement can 
be built only from the bipartite entanglement between 
fermion pairs. Lastly, we showed that the entropy of a 
large Fermi gas is always proportional to its volume, in- 
dependently of fermion distance. It is only for a small 
number of fermions and small distances that the entropy 
is smaller than this number (volume). We believe that 
our work shows that multipartite entanglement in com- 
plex macroscopic systems can be studied and even fully 
understood with the existing techniques of quantum in- 
formation. We hope that this stimulates other studies 
in similar directions of solid state and condensed matter 
systems. 
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